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1. SCOPE 


The "As-built" Design Specifications for the Production Film 
Production Film Converter Gains and Biases Program (PFCGAB) 
are the subject of this document. 

PFCGAB is a batch- input stand-alone program, programmed in 
the Univac Fortran V language, for operation under the EXEC 
II operating system on the Univac 1108 computer. 

The program was developed for the Earth Observations Division 

(TF) of NASA- JSC. The program requestor and source of 

* 

requirements was Mr. Richard Nance (TF3) . 


2 . APPLICABLE DOCUMENTS 


The following documents, of issue shown, are pertinent to the 
authorization for development and requirements for the 
PFCGAB program. 

• Requirements for a Computer Program to Investigate PFC 
Gains and Biases, by Richard Nance (TF3) (No date) 

• Action Document 2640: changes to PFCGAB, orignated by 

R. L. Nance (TF3) (No date) 

• LEC Interdepartmental Communication 643-1447, dated 9/28/76 

Subject: Modifications to Gains and Biases Program 


TIRF 76-0081 


3. SYSTEM DESCRIPTION 


A computer program, the Production Film Converter Gains and 
Biases program (PPCGAB) , is described in following sections. 

3.1 HARDWARE DESCRIPTION 

PFCGAB is operational on the Univac 1108 computer, with the 
EXEC2 operating system. The program may be made completely 
compatible with operation on the Univac 1110 and the EXECS 
operating system, with a simple conversion of the EXEC2 PCF 
(Program Complex File) to an EXECS PF (Program File) as docu- 
mented in the IDSD Procedures Manual, Part 20. 

3.2 SOFTWARE DESCRIPTION 

The program consists of a main element, PFCGAB, and twelve 
subprograms: IDHEAD, COVER, PRSTAT, PRXHIS, PFCHIS, PRPCT, 

PRGAB, TAPERD, UNPK8, UNPACK, TAPLAB, and BYTRAN. 

PFCGAB and the first 10 subprograms named are programmed in the 
Univac Fortran V language. The last two subprograms, BYTRAN 
and TAPLAB, are programmed in the Univac assembly language 
SLEUTHII. 

3.2.1 SOFTWARE COMPONENT NO.l (PFCGAB) 

PFCGAB (Production Film Converter Gains and Biases) is a pro- 
gram developed for the purpose of computing alternative gains 
and biases to be used by the PFC (Production Film Converter) 
when processing LACIE (Large Area Crop Inventory Experiment) 
images. The program is designed to input LACIE images from 
tape, and will accept either a GSFC (Goddard Space Flight 
Center) "CCT” LACIE image tape or ERIPS (Earth Resources 
Information Processing System) "merge" tape containing LACIE 
image (s). The program presumes 4-channels per image (ERTS 
acquisition) . 


3. 2. 1.1 Linkages 

PFCGAB is the main element. PFCGAB calls the following sub- 
programs : 

(1) An assembly language subprogram, TAPLAB, to provide the 
label of the input tape(s); (2) TAPERD, a Fortran V sub- 
progreun which reads the header information and image data 
from the input LACIE image t&pe^s); (3) PRXHIS,a Fortran V 
subprogram which prints the histogram (frequency distribution) 
of pixels (radiance values) in the LACIE image, (4) PRSTAT, a 
Fortran V subprogram which prints 4-channel statistics com- 
puted over the 5x6 nautical mile sample segment area and the 
10x11 nautical mile search area, (4) PRPCT, a Fortran V 
subprogram which prints a 4-channel tabulation of percentages 
of pixels, (6) PFCHIS, a Fortran V subprogram which prints 
a 4-channel histogram (frequency distribution) of PFC grey- 
shade levels (1-16) versus the number and percentage of 
pixels in the 5x6 nm sample segment area which are assigned 
to each grey-shade level, derived using the search area 
(10x11 nm) gains and biases Al, Bl, and (7) PRGAB, a Fortran 
V svibprogram which prints the 29 gains and biases derived in 
the main program PFCGAB. 

3. 2. 1.2 Interfaces 

PFCGAB interfaces with subprograms TAPLAB and TAPERD (entry 
points TAPHDR, FLDINT, and LINERD) by means of calling 
arguments. PFCGAB interfaces with the five output sub- 
programs PRXHIS, PRSTAT, PRPCT, PFCHIS, PRGAB - by named common 
blocks; ALL, PFCH, PRCENT, FREQ, HDATA, STCOMP, STlOll, 

ST5x6, PGLABL. 


3 . 2 . 1 . 3 Inputs 

Both card and tape input are required by PFCGAB. 


Tape input is required to be assigned to Uni vac logical unit A 
(Fortran unit 1) . Tape input is required to contain LACIE 
images either as provided on the LACIE (Large Area Crop 
Inventory Experiment) image tapes sent to JSC (Johnson Space- 
craft Center) from GSPC (Goddard Space Plight Center) on 
"CCT" ("Computer compatible Tape") or as provided from an 
ERIPS (Earth Resources Image Processing System) "merge" out- 
put of LACIE images to tape. The format of the tapes may 
either be Universal format or LARSYS II format. The program 
will accept and process multi-file tapes. 

Card input is expected by the program. Card input is required 
to be in Fortran V "NAMELIST" format. The "NAMELICT" name 
for the input is "$CARDIN." 

The parameters to be card- input are described in Section 4.0, 
Program Operation. 

3. 2. 1.4 Outputs 

PPCGAB utilizes the printer for all output. Output consists 
of five tabulations. All tabulations are 4-channel. The 
output is as follows: 

(1) A frequency distribution of pixel (radiance value) versus 
number and percentage of 5x6 nm sample segment area pixels 
having each radiance value, (2) basic statistics for the 5x6 
nm sample segment area (minimum radiance value, maximum 
radiance value, radiance value range, arithmetic mean (P2) , 
standard deviation (S2), mean deviation, median, and mode. The 
10x11 nm search area statistics output are the arithmetic mean 
(PI) and standard deviation (SI) , along with the search area 
gain (Al) and second bias (Bl) . (3) A tabulation of the per- 

centage of 5x6 nm sample segment area pixels which exceed a 
radiance value of 60, 90, 60-sun angle corrected, 70-sun angle 


corrected, 80-sun angle corrected, and 90-sun angle corrected* 
(4) A frequency distribution of the PPC grey-shade levels 
(1-16) versus the number and percentage of 5x6 nm sample seg- 
ment area pixels assigned to each PFC grey-shade level, using 
the search area gain and bias (Al, Bl) . (5) A print-out of 

all 10x11 nm search area and 5x6 nm sample segment area gains 
and biases derived by PPCGAB. 

If more than 4 channels are processed (i.e,, ERIPS "merge" 
images), the output tabulations are by 4-channel increments- 
i.e., channels 1-4, 5-8, 9-12, etc. 

All output is labeled such that the output variableis are self- 
identified and self-explanatory. 

3. 2. 1.5 Storage Requirements 

The computer core-storage requirement for PPCGAB and all sub- 
programs is 5956 I-bank locations, 16728 D-bank locations, 
and 2647 locations for all 9 Named common blocks. Total 
core-storeage requirement is 25331 locations. 

3. 2. 1.6 Description 

PPCGAB presumes that the input images are LACIE images as 
provided to JSC from Goddard Space Flight Center, or as 
derived on ERIPS and output on an ERIPS "merge" tape. The 
program considers an "image" to contain 4-channel pixels, 
with 196 pixels per scan line, and 117 "scan lines per "image." 
A maximum of 32 channels will be processed by the program; 
i.e., up to 8 "images," each "image" being a 4-channel image. 
This is the method used to process the ERIPS "merge" tapes 
with multiple LACIE acquisitions. 



PFCGAB inputs LACIE images from tape, computes statistics 
(mean, standard deviation) for the 10x11 nautical mile search 
area and the 5x6 nautical mile sample segment area (mean, 
standard deviation, range, mean deviation, median, emd mode)^ 

computes 5 gains and biases using 10x11 nm search area statis- 
tics, and computes 23 gains and biases using 5x6 nm sample 
segment area statistics. 


The 10x11 nautical mile search area gain is Al, and the second 
bias is Bl. A1 is either extracted from the input tape header 
for Goddard Tapes, or input via namelist "$CARDIN” for ERIPS 
"merge" tapes. Bl is computed from the search area bias, BO, 
extracted from the Goddard tape header or input via ”$SCARDIN.” 
Bl is computed as follows: 


Bli 



i = channel no. 


The 10x11 search area mean, PI, is computed as follows: 


PI. 

1 


"Bli ^^®/a 1 ^ 


The 10x11 search area standard deviation, SI, is computed as 
follows : 

Sl^ = 42.666662.,...—-^ 

— ^^i 

The five gains, biases computed using 10x11 nm search area 
statistics (PI, SI) are: A9, B9, AlO, BiO, A24, B24, A26, B26, A27, 

B27. The 10x11 nm search area related gains and biases are 
computed as follows: 

A9, B9 

a 5^ = 256/(40Sl^) 

B9. = (PI. - 2S1.) 

■1b Ja ^ 



i » 1, No. Channels 


AlO.BlO 

AlO. » 256/(5 SIJ 

1 1 

BIO^ = - (Pl^ - 2.5 Sl^) 

A24 , B24 (truncated) 

A24^ = 8 

B24j^ = (Plj^-16)*8 
A242 =5 

B242 = - (Pl2»26)*5 

A243 = 5 

624^ = - (Pl3-32)*5 

A24^ = 10 

B24^ = (Pl^-13)*10 

A2 6 > B2 6 (Bands 1,2,4 dependent) 

P new = MEAN = Pl^ + PI 2 + 2 Pl^ 

3 ^ 

M = MAX [(PI, + 3 Sl^ - (-PI. H. 3S1^ + P^^„ 

for i = channels 1, 2, and 4 
where for is* 4, Pl^ = 2*P1^ 

XHI = P „, + M 

new 

XLO = P , - M 

new 

A26 = 256/(XHI-XLO) 

B26 = - XLO * A26 

A27,B27 (Bands 2, 3, 4 dependent) 

= MEAN„^„ = P1_ + PI- + 2 pi. 
new new 2 3 4 

3 


H - MM] (Pl^ + 3Sli - , (-Pl^ + 3 Slj + P„^„) ) 

for i = 2, 3j» 4 an<l when i * 4, Pl^ 2*P1^ 

XHI = P + M 
new 

XLO = P „ 

new - M 

A27 = 256/(XHI - XLO), B27 » - XLO * A27 


The 5x6 nm sample segment area statistics are computed as 
follows: 


for X^ = radiance value; channel i 

NX^ = Number of pixels having radiance value X^; 

N^Otal “ Total number of pixels 

XMIN. = X . = MIN [X.] = minimum radiance value, channel i 

1 min. 1 


XMAX. = X = MAX [X.] = maximum radiance value, channel i 
1 max^ 1 


XRANGE. = RANGE. = X„,„ - X„. 

1 1 max^ minj 


N„ 


XMED^ = MEDIAN^ = X^ , for the X^ at which N^ > . TO^ AL 


XMODE. = MODE. = X. , for the X. satisfying MAX[N„ ] 

1 11 1 A. 


P2 . = MEAN . = 
1 1 


X. 


N„ 


JL. 


•TOTAL 

XMNDEV^ = MEAN DEVIATION ^ / ^ i " Z^i ( 


N 


TOTAL 


S2^ = STANDARD DEVIATION^ =1 



-1 









The gains and biases computed using the 5x6 nm sample segment 
area meems and standard deviation (P2 and S2) or mean deviation 
are as follows: 

A2 ,B2 


A2. * 256/6 S2. 

X X 


B2^ = - (P2^ - 3 S2^) 


A3,B3 (truncated distribution) 


P3^ 


X, 


= MEAN 


N 


TOTAL^ 


S3 



- 


for i = 2, 3< X.< 120 
— X— 

for i = 3, 5£ X^£ 120 
for i = A, 2£ X^£ 61 


A3^ = 256/6 S3^, B3^ = (P3^ - 3 S3^) 


A4 yB4 (truncated distribution) 


P4^ = MEAN 


S4 . = std. 

X 






where; for i * 1, 
i = 2, 
i = 3, 
i = 4, 


10 < < 80 

10 < < 80 

10 < Xi < 80 

5 < X^ < 40 


A4. = 256/6 S4., B4. = - (P4. - 3 S4.) 

X 1 1 1 1 

A5 ,B5 (truncated distribution) 


P5, 

X 


MEAN = 2 ^ 

N, 


S5 . = std. deviation 


TOTAL. 

.(s ) 

\ ®*TOTAL, / 



where , for i 
i 
i 
i 


= 1, 20 < X. < 60 

- X — 

=2, 20 < X. < 60 

— X — 

= 3, 20 < X^ < 60 
-4, 7 < X^ £ 30 


A6,B6 


A6j^ = 256/4 S2. 


B6. = (P2. - 2 S2 
X ' X 


hi ,B7 (uses mean deviation) 


A7^ = 256/6D^ 


where = mean deviation 


B7. = - (P2. - 3D, 

1 i ■ 

= El X. - P 2 <( 


N, 


TOTAL 


A8 yB8 (3a truncation) 


A8. = 

X 


256 


X 


MAX.- X„„. 

X MIN^ 


B8. = 
X 


- X 


MIN. 


d 


TOTAL^ 


) 


) 


where for N 

1 N 


X, 


MIN, 


X for N 


.00135 N^otAL “ upper tail) 

.00135 „ (LT « lower tail) 

“total 


A14 >B14 (2.5a truncation) 
256 


A14^ 


B14. 


- X 


^^i - ""min. 


MIN, 


where Xj^^ = ^UT„ ^ ” .0065 
1 N 

^MIN. “ ^LT ^ ^TOTAL 

X N 


A15,B15 (2.0a truncation) 


A15. = 256 

1 


^maxT ^min. 

X X 


B15 . — “ X„ym 
X MIN^ 


where — Xy^ for N .0225 

i N 


^MIN. “ ^LT„ ^ " *°225 N^Q^^y 

X N 


A16 ,B16 (3a truncation) 

A16, = 256 

X y m- V ^ 

“^MAX . MIN . 

X X 


B16 


- (P2^ 


^P2. '• 3a^ 

X 


MAX^"^ 

^UTn 

for 

N 

= .00135 

^TOTAL 

MIN^~ 

^LT 

for 

N 

= .00135 

N 

“total 

P2, “ 

3a = 

X. at 

X 

Xyp such 

that N^ 


’TOTAL 


where X^ is mean radiance value, band i 


jy{o 


A17»B17 (2.5o truncation) 


A17. 

3 

where > 
) 
) 

A18,B18 

A18. 

where 3 

4 

A2l>B21 

A21, 

where : 
with 


256 


X, 


, * “ ^^^i” ^P2^*‘ 2.5o) 


MAX^ " ^MIN, 


« X„„ for N » .00655 
MAX. UT„ ‘^TOTAL 

1 N 

^MIN. ^ * .00655 N^oTAL 

1 N 

nr- =® X . at X.. „ such that 

P2^ “ 2.5a i u.iai- 

N?s t N., ^ .49345 

^i ^i TOTAL 

where X^ » mean radiance, band i 
(2.0o truncation) 


256 


L " ^MAX. " ^MIN. 

1 1 


, B18^ " ^ 


P2. 


2a' 


“ X„_ for N *= .0229 
MAX. UT„ "total 

1 N 

MIN^ ® ^LTj^ ^ ^TOTAL 


P2^ - 2o = X^ at Xj^^ such that -f > .4771 N^oTAL 

where X. » mean radiance value, band i 

»!• 


(3a truncation) 



NEW, 




’NEW, 


NEW, 


E K, 


N„ 


TOTAL, 


constrained to be within the range 


P2. - 3 S2. < X. < P2. + 3 S2. 


and * Total no. pixels in the range P2£-3 S2^ < 

<P2^ 4 3 S2^ 


A22 ,B22 (iterated 3 j truncation) 



with X^ constrained within the range 

^ ®NEWj^ i i ^NEWj^ ^ ®NEWj^ 


and = total no, pixels in the range “ 

2 1 

X . < P„-,rT 4- 3 S„„„ 

1 • NEWj^ NEWj^ 




< 


A2 3 , B2 3 (truncated, sun angle corrected) 

A23i = f I - - . B23^ = - - 3 

= std. deviation of truncated distribution 
= mean of truncated distribution 

where ; 

for i = 1, 4*CFACT < X^ £ 80*CFACT 

i = 2, 4*CFACT < X^ < 80*CFACT 

i = 3, 4*CFACT < X^ £ 80*CFACT 




and = total no, radiance values, channel i, within the 

^ r 2 mge constraint for channel i 
and CFACT = SIN 60® 

A11,B11 


^^^i “ 5 S2 


SIN SUN ANGLE 
256 


Bll^ =*(P2^ - 2.5 S2^) 


A12,B12 


A32 = 


B12^ --(P2^ - 2.5 D^) 


where D. 

X 


, E 


X “ P2 

i ij = mean deviation, channel i 


N 


TOTAL 


A13,B13 


M3. = 

1 4D. 


, B13^ = - (P2^ - 2.0 D^) 


where = mean deviation in channel x 


. _ S I X . - P2 . 


N, 


TOTAL 



A25,B25 


A25j^ » 8 , B25j^ 
A252 = 5 , ®252 
A 253 = 5 , B 253 
A25^ = 10 , B25^ 


(P2j^-16)*8 

(P22-26)*5 

(P23»32)*5 

(P2.-13)*10 

4 


A19,B19 (Bands 1, 2, 4 dependent) 
256 

A19 = 


^MAX " ^MIN 


B19 = - X, 


MIN 


where * MAX [X^^l 

^MIN 


for i = Channels 1 , 2 , and 4, 

with constraint that X^ at i = 4 = 2*X^ 


A2 8 , B2 8 (Bands 1, 2, 4 dependent) 
256 


A28 = 


^HI " ^LO 


B28 = - Xj^q*A28 


where : 


^HI “ ^NEW ^ 

^LO ^NEW ~ ^ 

P^,„ = NEW MEAN = ^^1 ^ ^ ^^^4 


. M = MAX 


[(P2. 


+ 3 S2. - P.,„„ 
1 NEW 


) , (-P2. + 3 S2. + 

» X X 


NE 


for i = channels 1, 2 , and 4 


A20 ^B20 (Bands 2,3, 4 dependent) 
256 

^MAX “ ^MIN 


^MIN 




A20 


B20 


where = MAX [X.] 

haa 1 

=SlIN ' 

for i = channels 2 , 3, 4 and constraint X 
i = 4 


A29 /B29 (Bands 2 , 
256 


A29 = 


^HI " ^LO 


where Xjjj - 


^LO ” ^NEW 


^NEW = 


3, 4 dependent) 

, B29 = -Xj^q*A29 

+ M 
- M 

p2, + P2- + P2 
mean = — = »— = 


4 


M = MAX 


[( 


P2i+3S2.-P^^„ 


)'(- 


P2j^+3S2^+PNEW 


)] 


for i = channels 2 , 3 , 4 


PFCGAB also computes and outputs the frequency distribution of 
PFC grey-shade levels for the input image radiance values, using 
the 10x11 nm search area gain, Al, and second bias, Bl. The PFC 
grey-shade levels are computed such that there is equal pro- 
bability of image radiance Values being assigned to the 16 PFC 
grey levels, as follows: 

y. = Al. * (X.+Bl.) = gain, and bias corrected radiance value, 

^ ^ channel i 

^MAX. = MAX [Y. 1 


^MIN^ 


MIN [Y^] 


AZ 


^ = grey-level increment = — ^ 


16.0 


MAXj. - 


Ny = Number of Y-values = - Y 


‘MIN 


\lN^ 

+ 1.0 


LEVEL = J*AZ+1, J = l,2,...,Ny 

with constraint that 1(^ LEVEL) £ 16, integers. 
FREQUENCYj = J^^x ^ J . ' number of pixels at PFC 

level J = the nuni^er of pixels mapped into Y^ and Yj = 1*AZ 
< Yj < Yj + 1*AZ 

A sun angle may be input to PFCGAB from either the input image 
tape header (for GSFC tapes) or by NAMELIST $CARDIN card 
input. If available, a sun angle correction factor is com- 
puted as follows: 


CE’ACT = sun angle correction factor = 


SIN (60°) 

SIN (Sun angle) 


If the sun-angle is available the sun angle correction is 
applied to the 5x6 nm sample segment area statistics for 
output of sun angle-corrected statistics computed as follows: 

CXMIN. = CFACT * XMIN . = corrected X^^m 
1 1 MIN 

CXMAX. = CFACT * XMAX. = corrected X„.„ 

1 1 MAX 

CHANGE^ = CFACT * XRANGE^ = corrected range 


CMEAN^ 

= CFACT 

* 

P2^ « corrected mean 


CSTDEV^ 

» CFACT 

* 

S2^ * corrected standard 

deviation 

CMNDEV. 

1 

= CFACT 

* 

XMNDEV. * corrected mean 
1 

deviation 

CMODE , 
1 

CFACT 

* 

XMODE^ = corrected mode 


CMED^ 

* CFACT 

* 

XMED^ = corrected median 


for i 

s= band 1 

t 

2, 3, 4 


3. 2. 1.7 

Flowcharts 



The detail flowchart for PRCGAB is contained in Appendix A of 
this document. 

3. 2. 1.8 Program Listings 

The listing for every element in the PFCGAB program is con- 
tained in Appendix B of this document. 


4 . OPERATION 


Program PFCGAB is executed via a batch- run deck submitted to 
the NASA-JSC central computing complex in Building 12. The 
run deck must have as the innitial card in the deck the Univac 
"Run" card ("7/8 RUN ...") as documented in the Institutional 
Data Systems Division (IDSD) Procedures Manual, Part 19. The 
necessary Univac system control cards to assign the input data 
tape and the program tape to the run are illustrated in the 
"Sample Deck" illustration, section 4.5. Execution of the 
program, after the program tape contents have been entered 
into the Program Complex File area of the drum, is by means 
of the Univac Control card, "7/8 XQT PFCGAB." Following the 
"7/8 XQT PFCGAB" control card one or more Fortran V NAMELIST 
inputs are expected. The NAMELIST Name for input is "$CARDIN, 
with "$" always punched in card column 2 of the card. 

The final card in the batch-run deck, following all NAMELIST 
"$CARDIN. . . . $END" sets of variable values input to the pro- 
gram, should be the Univac End-of-File control card, "7/8 EOF. 
The EOF card is the program's signal to stop attempting to 
read "$CARDIN. . . $END" inputs, and terminate execution. 

The specific variables which may be card-input to the program 
via NAMELIST $CARDIN are described below. The user is 
cautioned that NAMELIST input demands that the variables be 
input by "Name = value," and that the Name must be spelled 
accurately in the input. 

The two required tapes to be input, and the program's printed 
output are also described below. 



4 . 1 CARD INPUT 


Input is by MEANS of NAMELIST "$CARDIN," i,e., ”$CARDIN 
VARl = XX, VAR2 = Y, VAR3 = Z, VAR4 (1) « XYZ.W, $END " 


The variables which may be input, ^ NAME , are: 


Variable 

GSFCTP 


REGTP 


Description 

(Integer) Flag to indicate that the input 
tape is a Goddard-generated tape- 
GSFCTP > 1, input tape ^ a Goddard tape 
GSFCTP < 0, input tape is not a Goddard 
tape 

(Integer) Flag to indicate that the input 
tape is one produced from segment- 
registration studies- 

REGTP > 1, input tape i£ a registration 
tape. 

REGTP < 0, input tape is not a "registra- 
ion” tape (NOTE; The capability to 
process registration tapes is not in the 

program at this time) 


MERGTP (Integer) Flag to indicate that the input 

Tape is a "merged-image” Tape from ERIPS- 
MERGTP > 1, input tape a merged-image 
tape . 

MERGTP < 0 , input tape is not a merged- 
image tape. 

NFILES (Integer) number of files to process, 

on the input tape (default = 1). 

FSTART (Integer) FIRST FILE NUMBER, at which 

processing is to begin on the input tape 
(default = 1) 

Processing of the input Tape begins at 


V x\ 


Variable 


NUMCH 


EL 


A 


B 


Description 

FSTART, and NFILES files will be sequen- 
tially processed on the Tape. 

(Integer) The total number of channels on 
the Tape-input image; should be an integral 
multiple of 4 channels (i.e.; 4, B, 12, 

16, 20, ..., 32) 

(default set = 4 channels) 

(Integer) Sun angle. If the Sun angle is 
input, the input sun angle will be used 
in computing sun- angle corrections to the 
5x6 area statistics. If the sun angle is 
not input, and is not available from the 
input tape header record, no sun-angle 
corrections will be attempted by the 
program. 

(default set = 0) 

(floating point) dimensioned = 4 - the 
10x11 search area gain for 4 channels 
(Al) , may be input as: ''A(l) = XX. X, 

A(2) = XX. X, A(3) = XX. X, A(4) = XX. X" 
(default set = 0.0) 

(floating point) dimensioned = 4 - the 
10x11 search area biases (Bl) for 4 
channels may be input as: ”B(1) - XX. X, 

B(2) = XX. X, B(3) = XX. X B(4) = XX. X" 
(default set = 0.0) 

(NOTE: If Al and Bl are input via A and 

B, the input gains and biases will be 
used, instead of the Tape-supplied Al, Bl, 
if available. If Al, Bl are 



Variable 


Description 

are neither input nor available 
from the tape, the gains and biases 
computed from Al, B1 are bypassed 
in the program, and will not be 
present in the printed output) 


The default (progiam-set) values assigned to the input 
variables are; 


GSFCTP = 1 
REGTP = 0 
MERGTP = 0 
FSTART = 1 
NFILES = 1 
MUMCH = 4 


(i.e., Goddard tape is the default) 


(file no. 1 is the default) 

(the number of files = 1 is the default) 
(total No. of channels on the tape = 4 is the 
default) 


A(l) , A(2) , A(3) , A4 = 0.0 
B(l) , B(2) , B(3) , B(4) =0.0 


Any variable not specifically input, ^ N7VME and value , will 
go into execution with the program-set default value. An 
"empty" NAMELIST input ("$CARDIN $END") will result 

in program execution with all input varjtsbles set to the 
default values. 

The program expects and will execute multiple "$CARDIN... 
$END" inputs. There is requirement to order multiple 
"$CARDIN $END" namelist inputs with respect to type of tape 
being processed, or the files to be processed. 


The program expects to be terminated by an "End-of-file” card 
("7/8 EOF"), following the last "$CARDIN. . . $END" NAMELIST 
input. 


4.2 TAPE INPUT 

The input Tape must be assigned to (logical) Univac unit "A” 
(Fortran unit 1). The input tape is presumed to be in the 
"Univeral" (or LARSYS2) tape format, one or more files, and 
four or more channels per input image up to a maximum of 32 
channels. The program internally processes "images" with 4- 
ciannels per "image" pixel, 196 pixels per scan line and 117 
■' c,n lines. The input tape is also presumed to be a 9 -TRACK 
tape, on the sample deck set-up shown in a succeeding section. 

The program (PCF) tape may be assigned to any unit other than 
"A" ("B" - "Z"). The program is entered into the computer 
from the program tape by appropriate use of the EXEC II system 
"CUR" ("Complex utility Routines") commands. See the sample 
deck set-up for the format and sequence of the control cards 
and commands to bring in the program from tape. 

The current Bldg. 12 library tape number for the program tape 
is X01400. 

4.3 OUTPUT 

All output is printer output, as follows; 

The first page contains various self-identified quantities 
from the "universal" tape header record, which are pertinent 
to the imagery in the file being processed. 

The second and succeeding pages of output are generated by 
program PFCGAB. Each page of output of a given type is 
headed by: 







File No. , of current input tape file being processed 
LACIE Sample Segment No., of current input tape file being 
processed. 

Sample Segment I.D . (LANDSAT Orbit) of the segment 
Latitude, Long.' >:ude , of saunple segment center 
Sun elevation . 

Sun Azimuth , at the time the imagery was LANDSAT-acguired 
(in degrees) 

All these quantities are obtained from the tape header record 
associated with the input tape file being processed. 

The first page of program PFCGAB output contains the un- 
corrected and (sun elevation) corrected statistics for the 
5x6 nm sample segment, and the statistics and bais and gain 
from the 10x11 nm search area, on a per-charmel basis. 

The second page of program PFCGAB output contains a tabulation 
of the percentage of pixels in the 5x6 nm sample segment, 
which exceed certain pre-selected pixel values, on a per- 
channel basis. 

The third (and probably fourth) pages of PFCGAB output are a 
frequency distribution histogram of the input 5x6 sample 
segment pixels, on a per-channel basis. 

The page following the 5x6 sample segment histogram is a 
histogram of the frequency distribution of the Production 
Film Converter (PFC) gray-shade levels for 5x6 nm sample 
segment, after application of the Goddard (or input) bias 
and scale factors (A1,B1) to the input imagery. 



Following the PPG Frequency Distribution, the PFCGAB program 
outputs the gain, bias, mean, and standard deviation for the 
10x11 nm search area (Al, Bl, PI, Sl)ythe 5x6 NM sample 
segment area (A2, B2, P2 , S2) and the program-computed gain, 
bias, mean, standard deviation for 3 pre-selected truncated 
distributions of the 5x6 NM area, the bias and scaling factor 
for "modified" (4 std dev.) standard deviation, the "average" 
(256/6X mean deviation) scaling and bias ( * sample mean - 3 
mean deviations) , and the "percent" scale and bias factor 
computed from the sample segment distribution within .00135 
tail limits. All computations shown on a per-channel basis. 
Succeeding pages of output contain the self-identified gains 
and biases computed by the program using various combinations 
of multiples of the standard deviation. See section 3. 2. 1.6, 
Description, for the method of computing the output gains and 
biases. 

4.4 PROGRAM USE RESTRICTIONS 

Although input is set up for three types of input tapes, the 
program is currently restricted to input and processing of the 
Goddard or ERIPS "merge" tapes only. Processing of the 
registration tapes requires additional specialized card-input 
variables to compensate for the lack of usable header infor- 
mation on the "registration" tapes. 

Also, the user is cautioned that NAMELIST input variables retain 
their values until specifically changed by input of the 
variable, by NAME, on a "$CARDIN. . . $END" namelist input. 

4.5 SAMPLE DECK SET-UP 

Card Column 1 

1 

7/8)zJRUN)zJ* • • 

7/8|rfSCH)zJ)zJ9T = 1, TT=1 


7/ 8 AN|rfASG)rf A=XXXXX, 

7/8bASG|i<B*X01400 
7/8)^XQT)rfCUR 
TRW A,B 
IN )6)if B 
TRI B 
TOC 

7/8)rf XQT PPCGAB 
Card Column 2 

r 

$CARDIN VAR. NAME1=X, VAR. NAME=XX, ... $END 
7/8)rfE0F 

4.6 PROGRAM EXECUTION CHARACTERISTICS 

Typical run-time, for a multi-file execution, is on the order 
of 1 minute per file. The amount (No. of pages) of print-out 
to be expected is estimated as 4 pages + no. files x 10 pages 
per file. 

The program executes on the UNIVAC 1108, EXEC2 operating 
system, in batch mode. The program run deck, user-supplied 
input tape, and the required form 588, "Instructions for 
Central Computer Complex Runs , " are submitted to the dispatch 
desk in building 12. The tape, if properly labeled by the 
user, and printed output from the run will be returned to the 
location code supplied on the "7/8 RUN" card of the deck, and 
noted in the proper place on the form 588 input to the dispatch 
area with the run deck. 

It is advisable to note on the I m 588, in the "Programmer 
Comments" area, that a 9-track tape unit is required for the 
run, and also on the form 588 under "Input Tapes," where the 
input tape unit and number are shown, to add a notation beside 
the input tape that it is a 9-track tape. 



NOTE: The program is not restricted to run with a 9-track in- 

put data tape. If LACIE imagery is to be input on a 7-track 
tape, the only required change is in the run deck, to use the 
correct Univac "7/8 ASG" control card options for the type 
of tape (9- track or 7- track) being input. 



APPENDIX A 


PFCGAB DETAILED FLOWCHART 








LAST 

FILE? 


READ PILE HEADER 
RECORD [TAPHDR] 

TO OBTAIN; 

SUN ELEVATION=SUNEL 

SUN AZIMUTH=SUNAZ 

BIAS, B (I), 1=IA CHANNELS 

GAIN, A°(I), 1^14 CHANNELS 

SITE I. D.=ISITE 

ERTS (LANDSAT) I.D. = IDERTS 

TIME OF IMAGE ACQUISITION; DAY, 

HOURS, MINUTE, TENS OF SECONDS = 
MISDAY, MISHR, MISMIN, MISEC 
SITE LOCATION (GEODETIC COORDINATES) ; 
DEGREES, MINUTE LATITUDE; LAT DEG, 

LAT MINi EAST/WEST; TEW 

DEGREES, MINUTES LONGITUDE; LON DEG, 

LON MIN; NORTH/SOUTH ; NORSOU 


SUN ANGLE (EL) 
INPUT ON CARDS 


FSUNEL=EL 


^ ^ .^NTYPE=1?^ 
(x.s. , GSFC 

Ifsunel^uneli-. 

compute X 

SUN ANGLE 1 
CORRECTION FACTOR ^ 
CFACT = SIN 60 “/SIN 
PSUNEL 


3 / 





SET FLAG [ELFLGl 

INDICATING PRESENCE OR ABSENCE 

OF SUN ANGLE 


COMPUTE TOTAL NO. OF 
IMAGES TO BE PROCESSED 
= NUMCH/4 


SET FLAG [AlBl] 

INDICATING PRESENCE OR 
ABSENCE OF 10x11 SEARCH 
AREA GAINS AND BIASES 
(Al, Bl) - IF GAINS, 

BIASES ARE INPUT USE 
THE INPUT GAINS AND BIASES 
IGNORE GAINS AND BIASES IN 
TAPE HEADER RECORD. IF 
GAINS AND BIASES NOT 
INPUT AND NOT AVAILABLE 
FROM INPUT TAPE HEADER 
RECORD, FLAG SET = .FALSE. 


READ AND 

HISTOGRAM DATA 

FOR 1 IMAGE (4 CHANNELS) 

PLACE IN IXHIST SUCH 

THAT IXHIST ( I, J)=FREQ OF 

X FOR CHANNEL I 

AND J=X=1, . . .1,256 



s 



FOR EACH BAND, i, 
i = 1/4, 

COMPUTE FOR 10x11 
NAUTICAL MILE SEARCH 
AREA; 

ARITHMETIC MEAN, Pi: 




=-Bl^ + 128/Al^ 


Standard deviation, SI; 
SI. = 42.666667/Al. 
10x11 najutical mile 
gains, biases: 

A9^ «= 256/4Sl^ 

B9^ =-(Pl^-2 Sl^) 

AlO^ = 256/5 Sl^ 

BIO^ =-{Pl^-2.5 Sl^) 

A24. = CA. 

1 1 


(CA = 8, 5, 


5, 10) 

B24^ = 

(CB = 16, 26, 32, 13) 


B24^ =-(Pl^-CB^)*CA^ 


(S 


> 4 " 

3 ^ 



ii 


SET Al, BX, Pi, SI, 

A9, B9, AlO, BIO = "BLANK” 

SET A24, B24, A26, B26, 

A27, B27 = BADUAL = 99999999.9 




COMPUTE GAIN, BIAS BASED 

ON MAXIMUM DATA VALUE FROM BANDS 1, 2, and 2* BAND 4 
AND MINIMUM DATA VALUE FROM BANDS 1, 2, and 2* BAND 4 
A19 = 256/ [MAX (BANDS 1, 2, 2*4) -MIN (BANDS 1, 2, 2*4)] 

B19 = -[MIN (BANDS 1, 2, 2*4) 


COMPUTE GAIN, BIAS BASED 

ON MAXIMUM DATA VALUE FROM BANDS 2, 3, and 2* BAND 4 
AND MINIMUM DATA VALUE FROM BANDS 2, 3, and 2 * BAND 4 
A20 » 256/ [MAX (BANDS 2, 3, 2*4) -MIN (BANDS 2, 3, 2*4)] 

B20 * - [MIN(BANDS 2,3, 2*4)] 


COMPUTE GAINS, BAISES BASED ON MEANS AND STANDARD 

DEVIATIONS (S2)„„„ DERIVED FROM DISTRIBUTION AFTER ELIMINATION 
OF POINTS MORE THAN 3*S2 FROM P2 ; 

A21 = 256/6 S2jjej^^ 

^ = -(P2j^g^^- 3 S2j,g^^^) 

A22 = 256/6 S2^^^ 

2 

B22 = - 3 S2jjg^j^) 


COMPUTE GAIN, BIASES FOR 10x11 NAUTICAL MILE SEARCH AREA, 
BASED ON SEARCH AREA MEAN, PI, AND STANDARD DEVIATION, SI, IN 
BANDS 1, 2, AND 4; 

A26 = 256/(H-L) 

B26 = - (L*A26) 


Where: P014 =(P1]^ 2*Pl^)/3 


M14 = MAX j (Pl^ + 3 Sl^-P014) , (-Pl^ + 3 

i = BAND 1, 2, 2*4 

H = P014+M14 

L = P014-M14 





m 








m 




J = J+1 


I SUMX 
SUMXSQ 
FNUM 

SUMSQ 

XMEAN 

VAR 

STDDEV 


MAX 


IJ 


IXHIST(J,K)*K 

K=MIN, 


IJ 


MAX 


E l J f 

IXHIST(J,K)*K" 

MAX 


IJ 

NX = / . IXHIST(J,K) 

K=MIN_ _ 

lu 

NO. PIXELS, FOR CHANNEL J 

SUMX* SUMX 

SUMX/FNUM 

VARIANCE = (SUMXSQ-SUMSQ/FNUM)/(FNUM-1.0) 
STANDARD DEVIATION = n/VAR 





P3(J) - XMEAN 
S3(J) » STDDEV 
A3(J) - 256/[6*S3(J)] 
B3(J) * -IP3(J)-3*S3(J)1 


P4(J) » XMEAN 

S4(J) » STDDEV 

A4(J) * 256/l6*S4(J)l 

B4(J) * - [P4(J) - 3*S4(J)] 


(truncate 
^ang 


P5(J) * XMEAN 

S5(J) = STDDEV 

A5(J) * 256/[6*S5(J)] 

B5(J) = - [P5(J) - 3*S5(J)1 


P23(J) = XMEAN 
IT '|S23(J) = STDDEV 

A23(J) = 256/l6*S23(J) ] 
■IeLFLgS Ib23(J) = - [P23(J) - 3*S23(J)] 


P23 (J) 

= 99999999.9 

(BADVAL) 

S23(J) 

= 99999999.9 

(BADVAL) 

A23 (J) 

= 99999999.9 

(BADVAL) 

B23 (J) 

= 99999999.9 

(BADVAL) 


COMPUTE 3o, 2.5a, 2.0a GAINS AND BIASES 
(A8, B8, A14, B14, A15, B15, A16 , B16 , 
A17, B17, A18, B18) : 


NUMX^ = TOTAL NO. PIXELS IN CHANNEL i 


— GSCALEj^ / k — \ , . . . ,5 

= .00135, .00655, .0225, .00135, .00655, 
.0229 

= BSCALEj^, k = 1, . . . , 6 

= .00135, .00655, .0225, .49865 
.49345, .4771 





FOR COMPUTATION OP GAINS FOR CHANNEL i ; 
MIN. = N (PIXEL VALUE IN CHANNEL i) 

= M-1 FOR M SUCH THAT 


IXHIST(J). > GSC.*NUMX. 

X K X 


MAX.= N (PIXEL VALUE IN CHANNEL i) 

= M-1 FOR M SUCH THAT 
_2_56 

2 IXHIST(J)^ > GSCj^*NUMX^ 


for k = 1 , 2 , . . . , 6 



FOR COMPUTATION OF BIASES FOR CHANNEL i; 
FOR k = 1, 2 , 3 : 

MIN. = M-1 FOR M SUCH THAT 
^k M 

;rXHIST(J)^ > BSCj^*NUMX^ 


FOR k = 4, 5, 6; 


MIN. = M-1 FOR M SUCH THAT 

^k P2j^ 


IXHIST(J)j^ ^ BSCj^*NUMX^ 






A17. = 256/ (MAX. - MIN. ) 


= 2.5a GAIN, CHANNEL i 

A18. = 256/ (MAX. - MIN. ) 

^ ^6 ^6 

3.0,2.5,2.0a BIASES;'^ GAIN, CHANNEL i 

B8. = “ MIN. 

B14. = “ MIN. 

2 

B15. = - MIN, 

3 

B16. = - MIN. 

X 

B17. = -MIN. 

1 


B 18 ^ = - MIN^ 




CXMAX.= ~ CORRECTED MAXIMUM PIXEL VALUE, 

1 sinTeiT max. cahnnel i 








GRANGE. = j * (X„. „ - ) = CORRECTED RANGE, 

1 SIN (EL) MAX^ MIN^ CHANNEL i 

CMEAN^ = P2^ = CORRECTED MEAN 

CSTDEV^ = fiiriiT* “ CORRECTED STANDARD DEVIATION 


CMNDEV^ 


XMNDEV^ = CORRECTED MEAN DEVIATION 


CMODE . 


SIN (60)* XMODE. = CORRECTED MODE 

sinTelT ^ 


CMED. 


= * XMED^ = CORRECTED MEDIAN 


FREQUENCY DISTRIBUTION OF PFC DATA VALUES, y; 

FOR lY = Yi = Al^ * (X^ + Blj^) +0.5, 

lYHIST (i, lY, 1) = IXHIST (i, IX) 

= NO. OF PIXELS, CHANNEL i 

PFC GREY-LEVEL ASSIGNMENT (EQUAL PROBABILITY OF Y-VALUES 
ASSIGNED TO PFC GREY- LEVEL) : 


LVLPPC 


PFC GREY -LEVEL, 1< LVLPFC<16 


AZ = GREY-LEVEL INCREMENT 


=. 16.0 

Y -Y 
MAX MIN 


NO. GREY LE VELS 
Y RANGE 


FOR J2 = 1, 2, 3, . . .Ny 

AND Ny = Yj^ - +1.0 NO. OF Y-VALUES, 

IZ = J2*AZ+1.0 





WITH CONSTRAINTS: 

IF IZ>16, IZ=16 
lYHIST (i, lY, 2) = IZ 

= PPG GREY-LEVEL FOR Y = lY 


FREQUENCY DISTRIBUTION OP PPC GREY- LEVELS: 

HISPFC (i, IZ, 1) = HISPFC(i, IZ, 1)+ lYiilST (i, J2, 1) 

= NO. OF Y- VALUES MAPPED INTO PFC 
GREY- LEVEL IZ, 1<IZ^16 


HISPFC (i,IZ, 2) = HISPFC (i,IZ,l) 

y TOTAL 


* 100.0 


% Y- VALUES IN PFC GREY LEVEL IZ 



(REGISTRATION) =‘2 


IMAG 

sTYPE 


=1 (GODDARD TAPE) 

=3 (ERIPS ''MERGE” TAPE) 


NOT 

AVAILABLE- 
GET INPUT 
FOR NEXT 
CASE 


PRINT PAGE HEADING: 


PAGE HEADING: 


PAGHD2, WITH VARIABLES PAGHDl, WITH VARIABLES 
NFPROC, TLABEL, NIMAGE, TLABEL, LATDEG, LATMIN, 


PSUNEL 


NORSOU, NFPROC, ISITE, 
IDERTS, MISDAY, MISHR, 
MISMIN, MISSEC, LONDEG, 
LONMIN, lEW, PSUNEL, 
SUNAZ 


PRINT FREQUENCY DISTRIBUTION 
FOR 5x6 NM SAMPLE SEGMENT AREA: 

FOR EACH CHANNEL, 

PIXEL VALUE VS NO. OF PIXELS OF THAT 
VALUE, AND PERCENTAGE OF AREA HAVING 
PIXEL VALUE 

THAT 

PRINT STATISTICS BY CHANNEL, FOR 5x6 NM 
SAMPLE SEGMENT AREA, AND 
10x11 NM SEARCH AREA, WITH AND WITHOUT 
SUNANGLE CORRECTIONS: 

5x6 AREA: 

^MiN' ^MAx' Grange' 

Xa(S2), XMa (MEAN DEVIATION), 


^MEDIAN' ^MODE' 


10x11 AREA: 


Al(GAIN), Bl (SECOND BIAS) 

PI (MEAN), SI (STD, DEVIATION) 



A- <E 




10x11 NM 
SEARCH AREA 
GAIN (Al) AND SECOl^ 
BIAS (Bl) 

J^VAILABLE? 


^BS 













APPENDIX B 


PFCGAB PROGRAM LISTINGS 
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